Theory of quantum impurities in spin liquids 
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We describe spin correlations in the vicinity of a generalized impurity in a wide class of fraction- 
alized spin liquid states. We argue that the primary characterization of the impurity is its electric 
charge under the gauge field describing singlet excitations in the spin liquid. We focus on two gapless 
U(l) spin liquids described by 2+1 dimensional conformal field theories (CFT): the staggered fiux 
(sF) spin liquid, and the deconfined critical point between the Neel and valence-bond-solid (VBS) 
states. In these cases, the electric charge is argued to be an exactly marginal perturbation of the 
CFT. Consequently, the impurity susceptibility has a 1/T temperature dependence, with an anoma- 
lous Curie constant which is a universal number associated with the CFT. One unexpected feature 
of the CFT of the sF state is that an applied magnetic field does not induce any staggered spin 
polarization in the vicinity of the impurity (while such a staggered magnetization is present for the 
Neel- VBS case). These results differ significantly from earlier theories of vacancies in the sF state, 
and we explicitly demonstrate how our gauge theory corrects these works. We discuss implications 
of our results for the cuprate superconductors, organic Mott insulators, and graphene. 

PACS numbers; 75.30.Hx,75.40.Cx,76.60.Cq,74.25.Ha 



I. INTRODUCTION 

The response of a strongly interacting electronic sys- 
tem to impurities has long been a fruitful way of experi- 
mentally and theoretically elucidating the subtle correla- 
tions in its many-body ground state wavefunction. The 
most prominent example is the Kondo effect, which de- 
scribes the interplay between a variety of impurities with 
a spin and/or 'flavor' degree of freedom and a system of 
free fermions with either a finit o"'^'^'^ or vanishing^ density 
of states at the Fermi energy. 

More recently, the impurity responses of a va- 
riety of 'non-Fermi-liquid' bulk states have been 
studied,'^"i7iVd"dbi2.i3,i4,i5,i6,i7,i8 The S = 1/2 anti- 

ferromagnetic spin chain generically has a critical ground 
state, and displays interesting universal characteristics in 
its response to impurities or boundaries^. Universality 
was also found in the general theorjiSii^ of impurities 
in 'dimerized' quantum antiferromagnets in spatial di- 
mensions d > 2 near a quantum critical point between a 
Neel state and a confining spin gap state. Such 'dimer- 
ized' antiferromagnets have an even number of = 1/2 
spins per unit cell, and consequently their bulk quantum 
criticality is described within the conventional Landau- 
Ginzburg- Wilson (LGW) framework of a fluctuating Neel 
order parameter. ^^■^''•^^ Away from the impurity, such 
systems only have excitations which carry integer spin. 

It is the purpose of this paper to extend the above 
theorji£i24i& to fractionalized 'spin liquid' states in spa- 
tial dimensions d > 2 with neutral S — 1/2 excita- 
tions ('spinous') in the bulk. Such spinon excitations 
carry gauge charges associated with an 'emergent' gauge 
force (distinct from the electromagnetic forces) , typically 
with the gauge group^^'^'^ Z2 or U(l), and we will ar- 
gue shortly that such gauge forces play a key role in 
the response of spin liquid states to impurities. Earlier 
analyses^^iSS*2SiSL2& of the influence of impurities in the 



U(l) 'staggered-flux' spin liquid ignored the crucial gauge 
forces; we will comment in detail on the relationship of 
our results to these works in Section Fill Bl An analysis 
of impurities in a Z2 spin liquid state was presented re- 
cently by Florens et al.f^ but for a particular situation 
in which a spin moment was strongly localized on an im- 
purity and gauge forces could be safely neglected; we will 
comment further on their work in Section [IV Bl 

There are a number of experimental motivations for 
our work. A large number of experiments have studied 
Zn and Ni impurities in the cupratest22i22i^iiS and much 
useful information has been obtained on the spatial and 
temperature dependence of the induced moments around 
the impurity. It would clearly be useful to compare these 
results with the corresponding predictions for different 
spin liquid states, and for states proximate to quantum 
critical points. We will show here that there are signifi- 
cant differences in the experimental signatures of the dif- 
ferent candidates, and this should eventually allow clear 
discrimination by a comparison to experimental results. 
A second motivation comes from a recent nuclear mag- 
netic resonance (NMR) studji^^ of the 5 = 1/2 trian- 
gular lattice organic Mott insulator k-(ET)2Cu2(CN)3, 
which possibly has a non-magnetic, spin singlet ground 
state. The NMR signal shows significant inhomogeneous 
broadening, indicative of local fields nucleated around 
impurities. Our theoretical predictions here for Knight 
shift around impurities should also assist here in selecting 
among the candidate ground states. 

An important observations is that in situations with 
deconfinement in the bulk, the bulk spinous are readily 
available to screen any moments associated with an im- 
purity atom (as has also been noted by Florens et al}^). 
Moreover, for a non-magnetic impurity (such as Zn on a 
Cu site), there is no a priori reason for the impurity to 
acquire a strongly localized moment. Consequently, it is 
very useful to consider the case where the impurity has 
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local net spin S = 0. Naively, such a situation might 
seem quite uninteresting, as there is then no local spin 
degree of freedom which can interact non-trivially with 
the excitations of the spin liquid. Indeed, in a Fermi 
liquid, a non-magnetic impurity has little effect, apart 
from a local renormalization of Fermi liquid parameters, 
and there is no Kondo physics. However, in spin liquids 
the impurity can carry an electric gauge charge Q, and 
the primary purpose of this paper will be to demonstrate 
that a Q ^ 0, S — impurity displays rich and universal 
physics. 

We will primarily consider U(l) spin liquids here: then 
an important dynamical degree of freedom is a U(l) 
gauge field Afj,, where extends over the d + 1 space- 
time directions, including the imaginary time direction, 
T. Our considerations here can also be easily extended 
to Z2 spin liquids, and this will be described later in 
Section liVBI We normalize the gauge field such that 
the spinous have electric charges ±1. In the popular U(l) 
gauge theories of antiferromagnets on the square lattice 
(which we will describe more specifically below), a va- 
cancy will carry a gauge charge Q = ±1. Thus, a Zn 
impurity on the Cu square lattice site has Q = ±1. This 
can be understood by thinking of the impurity as a lo- 
calized 'holon' in the doped antiferromagnet, which also 
carries such gauge charges 

We will consider theories here with actions of the struc- 
ture 

S = Sb + Simp, (1-1) 

where Sb is the bulk action of the spin liquid in the ab- 
sence of any impurity, and Simp represents the pertur- 
bation due to an impurity which we assume is localized 
near the origin of spatial coordinates, x — 0. We will 
argue that the dominant term in Simp is the coupling of 
the impurity to the U(l) gauge field: 

Simp = iQ J dTAr{x^O,T). (1.2) 

We will demonstrate that additional terms in the im- 
purity action are unimportant or 'irrelevant'. The iSimp 
above can be regarded as the remnant of the spin Berry 
phase that characterized the impurity in the previous 
theory2i24i& of dimerized antiferromagnets; the latter 
Berry phase for a spin S impurity was iS times the area 
enclosed by the path mapped on the unit sphere by the 
time history of the impurity spin. An explicit reduc- 
tion in the spinon formulation of the spin Berry phase to 
Eq. H1.2|l was presented in Ref. 's^. 

We will consider a variety of realizations of the bulk 
spin liquid Sb in this paper. Our primary results will be 
for U(l) algebraic spin liquids, ■^^ which are described by 
2-1-1 dimensional conformal field theories (CFT) and we 
specialize our presentation to these CFT cases in the re- 
mainder of this section. An algebraic spin liquid has gap- 
less spinon excitations which interact strongly with the 
gauge field. An explicit realization appears in the de- 
confined quantum critical pointii*2Si2i between Neel and 



valence bond solid (VBS) states, in which the spinous 
are relativistic bosons described by the CP^~^ field the- 
ory. Another is found in the 'staggered flux' (sF) phase 
of SU(Af) antiferromagnets, where the spinous are Dirac 
fermions'^^''^^''^^''^*. In all these cases, the algebraic spin 
liquid is described by a 2-fl dimensional conformal field 
theory, and our primary purpose here is to describe the 
boundary conformal field theory that appears in the pres- 
ence of 5i,np. 

Our central observation, forming the basis of our re- 
sults, is that 5iinp in Eq. H1.2|l is an exactly marginal 
perturbation to the bulk conformal field theory. This 
non-renormalization is a consequence of U(l) gauge in- 
variance, which holds both in the bulk and on the im- 
purity. We will verify this non-renormalization claim in 
a variety of perturbative analyses of the conformal field 
theory. The claim can also be viewed as a descendant 
of the non-renormalization of the spin Berry phase term, 
found in Ref. 

The exact marginality of iSimp has immediate conse- 
quences for the response of the system to a uniform ap- 
plied magnetic field H . The impurity susceptibility, ximp, 
defined as the change in the total bulk susceptibility due 
to the presence of the impurity, obeys 

Ximp = ^ (1-3) 

at finite temperature T above a conformal ground state; 
this can be extended by standard scaling forms (as in 
Ref. 9^ to proximate gapped or ordered phases, as we will 
describe in the body of the paper. We set h — ks — ^ and 
absorb a factor of the magneton, g/is , in the definition 
of the Zeeman field. With this, C is a dimensionless uni- 
versal number, dependent only upon the value of Q, and 
the universality class of the bulk conformal field theory. 

It is remarkable that the response of the impurity has a 
Curie-like T dependence, albeit with an anomalous Curie 
constant C (which is likely an irrational number). This 
anomalous Curie response appears even though there is 
no spin moment localized on the impurity. In contrast, 
the earlier results for the LGW quantum critical point 
presented in Ref. J, had an unscreened moment present 
and so a Curie response did not appear as remarkable. 
Here, it is due to the deformation of a continuum of bulk 
excitations by the impurity, and the 1/T power-law is a 
simple consequence of the fact that H and T both scale 
as an energy. Indeed any other external field, coupling to 
a total conserved charge, will also have a corresponding 
universal 1 /T susceptibility. 

A Curie-like response of an impurity in the staggered 
flux phase was also noted early on by Khalliulin and 
collaborators and othersiSSiSi However, in their mean- 
field analysis, they associated this response with a zero 
energy 'bound state', and hence argued that C = 1/4. As 
noted above, the actual interpretation is different: there 
is a critical continuum of excitations, and its collective 
boundary critical response has a Curie temperature de- 
pendence as a consequence of hyperscaling properties. 
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Consequently, C does not equal the Curie constant of a 
single spin, and is a non-trivial number which is almost 
certainly irrational. We will discuss the earlier work more 
explicitly in Section IlII Bl 

We will also consider the spatial dependence of the re- 
sponse to a uniform applied field, H, in the presence of an 
impurity, as that determines the Knight shift in NMR ex- 
periments. The uniform magnetization density induced 
by the applied field leads to a Knight shift HKu{x); at 
T above a conformal ground state this obeys (the scaling 
form is also as in Ref. 0): 

Ku{x) = ^-±^^^{xT/c) (1.4) 

where c is the spinon velocity in the bulk (we assume the 
bulk theory has dynamic critical exponent z = 1, and 
henceforth set c = 1), <!>„ is a universal function, and the 
Knight shift is normalized so that 

j (fxKu{x) = Ximp- (1-5) 

The function $«(?;) has a power-law singularity as y — > 0, 
with the exponent determined by a 'boundary scaling 
dimension': this will be described in the body of the 
paper for the various cases. 

In addition to the locally uniform Knight shift, an im- 
purity in the presence of a uniform applied field also in- 
duces a 'staggered' moment which typically oscillates at 
the wavevector associated with a proximate magnetically 
ordered state. This leads to a staggered Knight shift, 
HKs{x), which we will also consider here. Such a stag- 
gered Knight shift does appear for the deconfined criti- 
cal theory describing the Neel-VBS transition, and it has 
a spatial distribution associated with that of the Neel 
state. However, the response for the U(l) sF spin liq- 
uid is dramatically different. One of our primary results 
is that for the scaling limit theory of the U(l) sF spin 
liquid, an applied magnetic field in the presence of an 
impurity induces none of the many competing orders'^'' 
associated with the spin liquid. Thus there is no analog 
of the 'staggered' Knight shift. A subdominant induction 
of competing orders can arise upon including irrelevant 
operators associated with corrections to scaling; the pri- 
mary response, however, is just the induction of a fer- 
romagnetic moment, which has a slowly- varying, space- 
dependent envelope in the vicinity of the impurity spec- 
ified by Ku{x). 

Our conclusions above for the impurity response of 
the U(l) sF spin liquid differ from the earlier mean- 
field theories. They found an induced moment 
which had a strong oscillation between the two sublat- 
tices of the square lattice. We demonstrate here that 
this oscillation disappears in the continuum field theory 
which accounts for the gauge fluctuations. We are not 
aware of any reason why fluctuation corrections to the 
mean-field predictions should be considered small. 

The outline of the remainder of the paper is as follows. 
In Section^ we will consider the CP^~^ model field 



theory which describes the vicinity of the Neel-VBS tran- 
sition. The U(l) sF spin liquid will then be described in 
Section UTTI Other spin fiquids, not described by CFTs 
or by exactly marginal impurity perturbations, will be 
briefly discussed in Scction llVl Finally, experimental im- 
plications and a summary appear in Section IVI 

II. CP'^-^ MODEL 

The CP^^^ model describes the deconfined quantum 
critical point between the Neel and VBS states on the 
square latticeii^iSSiSi It is a field theory of complex scalars 
Za, a = 1...N, with global SU(A^) symmetry and a 
coupling to the U(l) gauge field A^. Near the transition 
at which global SU(7V) symmetry is broken we can work 
with the effective action 




where y = (r, x) is a spacetime coordinate, extends over 
D spacetime indices, D is the dimension of spacetime, 
related to the spatial dimensionality, d, by 

D = d+1 (2.2) 

We are interested in the phases of the field theory in 
Eq. H2.1() as a function of the tuning parameter s. For 
■3^0 there is a magnetically ordered phase with global 
SU(Af) symmetry broken, while for s :s> we have a 
spin-gap phase with full SU(Af) symmetry. We are es- 
pecially interested in the conformally-invariant critical 
point which separates these phases. We will begin by re- 
viewing the critical properties of the bulk action Sb alone 
in the following subsection, and describe impurity effects 
in the subsequent subsections. 

A. Bulk theory 

We will restrict our analysis here to the e expansion 
where 

e = 4-D. (2.3) 

This was carried out by Halperin et al.^^ some time 
ago, in a different physical context. This section will 
merely recast their results in our notation, using the field- 
theoretical formulation. 

As noted by Halperin et ai, a stable fixed point is ob- 
tained in the e expansion only for sufficiently large values 
of TV. However, it is expected that in the physical dimen- 
sion of e = 1, the fixed point may well be stable down 
to the needed values of N. Our purpose here is to un- 
derstand the basic features of the second-order critical 
point, and its response to impurities: so we will assume 
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that the value of N is large enough to ensure stability of 
the fixed point. This will delineate the essential scaling 
structure of the impurity response, but is not expected 
to be quantitatively accurate. 

It is also possible to analyze this model using the 
expansion, directly in D = 3. We choose not to present 
this here, because the results are very similar to the e- 
expansion, and the methods are closely related to those 
used for the 1/Nf expansion in Section UTTl 

The renormalization proceeds by defining rcnormalized 
fields and coupling constants u and / by 

— 7I/2 



Uo 



ZISd 



Sr 



f 



where 



Sn — 



r(D/2)(47r)^/2 



(2.4) 



(2.5) 



We work in the Lorentz gauge, in which the propa- 
gator is {Sp_i, — p^pv/p^)/p^. In this gauge the renormal- 
ization constants are 



1 



1 



{N + 4)7i 

El 
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(2.6) 



From this we can determine the anomalous dimension, 
rjz of the Za] note that this is gauge dependent 



•Hz 



= -3/, 



(2.7) 



where this is to be evaluated at the fixed point of the /3 
functions in Eq. (|2.8|l below. 
The (3 functions are 



Piu) 
/?(/) 



= -eu +(N + 4)u2 + Qf - Qfu 

dfi 



d^ 3 



(2.8) 



The renormalization of determines the critical ex- 
ponent v. This is associated with the renormalization 
constant 



Zo = 1 



{N + 4)u 



(2.9) 



In this gauge, there is no contribution to Z2 of order /. 
The critical exponent v is given by 



1 d Z2 

V dfi Zz 

= 2-{N + A)u + 3/, 



(2.10) 



where, again, this is to be evaluated at the fixed point of 
the P functions in Eq. (|2.8|) . 

Finally, we consider the scaling dimensions of gauge- 
invariant operators which characterize the observable 
spin correlations. The Neel order parameter is defined 
by 



(2.11) 



where is a iV x matrix which is a generator of 
SU(A^). We define its anomalous dimension -q by 



ba]^{D-2 + ri)/2 



(2.12) 



In the Lorentz gauge, (pa has no additional renormal- 
ization from A^^ fluctuations at leading order in e. The 
anomalous exponent rj is then given, to this order, by 



rj ^ D -2 + 2riz + 2u + 0{e'^) 



(2.13) 



Note that the value of rj is gauge-invariant, while that of 
rjz is not; the relationship H2.13|l holds only in the Lorentz 
gauge. 

We will also be interested in correlations of the mag- 
netization density Ma- This is defined by the response 
of the system to a uniform and staggered magnetic fields 
Hu and , under which the action is modified by 

x{(9, - iAr)z,, - KT^^z^} - H^Pa. (2.14) 
The magnetization density is given by 

SlnZ 



Ma = T 



5H„, 



(2.15) 



where Z is the partition function. Because Hu couples to 
a conserved 'charge', it scales as an energy, and therefore 



dim[Ma] = d. 



B. Impurity exponents 



(2.16) 



We now turn to an analysis of spin correlations in the 
vicinity of the impurity. The general method is very sim- 
ilar to that followed in Ref. UM 

We assume here that monopole tunneling events re- 
main irrelevant at the impurity at the quantum critical 
point, as they do in the bulkii2iS2i2i Such monopole tun- 
neling events are defined only in d = 2, and so their 
scaling dimensions are not easily estimated in the e ex- 
pansion. However, the monopoles are irrelevant at large 
N in the bulk because their action is proportional to TV, 
and the same reasoning applies also in the presence of 
the impurity. 

As in Refs. I^ITiL the staggered and uniform magneti- 
zations operators {(pa and Ma) acquire additional renor- 
malizations as they approach the impurity at x = 0. Let 
us denote the corrections by modifying H2.12|) to 



dim[Pa{x 0)] = dim[Pa] + 



imp 



(2.17) 



and similarly for Ma- We will compute these additional 
renormalizations to order below, and find that to this 
order A?^ = A*^ 

U1U.C1 ^iinp "imp- 

It is also useful to restate the above results in the no- 
tation of Refs. Alfk These we find that the excitations 
near the impurity are characterized by a spin operator 
Sa^ where a is the index of SU(Af) generators, and we 
denoted its scaling dimension by ?7'/2, or, equivalently, 



(5a(T)5b(0)) 



5ab 



(2.18) 



at the quantum critical point. Assuming that Sa is ob- 
tained as (j)a approaches the impurity, we obtain from 
and ifTTTIl that 



r/ = i?-2 + r; + 2Af^p. 



(2.19) 



These scaling relations imply*'' that at the quantum crit- 
ical point the temperature dependence of the NMR re- 
laxation rate is 



1/Ti - T-^+^' 



(2.20) 



near the impurity, while l/Ti ^ j'd-2+-q j^-^ ^y^q bulk. Fi- 
nally, we note that we can also express the above relations 
using the the operator product expansion 



lim 4>a{x,T) 



Sa{r) 



(2.21) 



It is also useful to consider time-dependent correlations 
of the magnetization density, Ma , as one approaches the 
impurity. The initial guess would be that these acquire 
the additional anomalous dimension A^^p. However, we 
will see in our computations below that the dominant 
correlations are instead given by 'mixing' between the 
4>a and Ma operators near the impurity. As was also 
found in Ref. |9|, the impurity Berry phase allows such 
mixing, and so the operator product expansion of Ma 
near the impurity has a dominant term given by the same 
operator Sa introduced above to describe the behavior 
of <j)a- This assertion is encapsulated in the operator 
product expansion 



lim Ma{x,T) 



Sa{r) 



(2.22) 



The relations H2.18|l . (|2.21() and (|2.22() are keys to our 
analysis, and will allow us to deduce the structure of 
spin correlations near the impurity. Their validity will 
be supported by a number of perturbative results we will 
obtain below. 



Now we turn to a determination of the exponents A? 



imp 



andAj^j^jp. These are most easily specified by determining 
the additional renormalization factor needed to cancel 
poles in e in correlators of (j)a and Ma as one approaches 
the impurity. The simplest correlator of (pa which allows 
us to achieve this is the vertex operator 

{(j)a{x = 0, w = 0)za(fc = 0, uj)zl{k = 0, Uj)) , (2.23) 




FIG. 1: Feynman diagrams which contribute to the impurity- 
dependent renormalization of vertices such as those in 
Eq. 12.231 . The full line is the Za propagator, the wavy line 
is the propagator, and the filled square is the source term 
in (Simp (also in Fig. I^Ib) . 



where the external Zq. legs will be truncated, and simi- 
larly for the Ma ■ Here the external frequency w is needed 
to control the infra-red singularities. 

The leading contribution to (|2.23|l , at order e is shown 
in Fig.^a) and it renormalized the vertex by the factor 



q 



1 



(27r)'^92(q2+^2 



(2.24) 



The same expression also applies to the renormalizations 
of the vertex associated with Ma- The q integral is easily 
evaluated, and no pole in e appears; this is fortunate, 
because the expression in (|2.24|l is purely imaginary. 

The determination of the O(e^) corrections involves 
a more difficult computation. There are a number of 
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FIG. 2: Vertices appearing in the calculation of response 
to applied field; Hu and Ha denote uniform and staggered 
magnetic fields, respectively. 
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FIG. 3: Feynman diagrams determining the leading order 
contribution to the free energy in presence of uniform (dia- 
grams (a)- (g)) and staggered (diagrams (h), (i)) magnetic 
fields. 



graphs which contribute, shown in Fig. ^ Fig. ^b)) is 
simply the square of Fig. ^a), and so contributes no pole 
in e. The value of Fig. ^c) (and its symmetry related 
partner) is 



1 



{2ttY (27r)'^ (72(^2 
no pole in e. 



-u;2)(p 



(2.25) 



The contribution of Fig. (and its symmetry related 
partner), for both the 4>a and Ma vertices is 



2Q 



d'^q d'^k 



1 



{2ttY {2nY q^(k^ +uj^){q + kY 
4^2/2 



(2.26) 



The diagrams Fig. ^e-j) are all purely imaginary, so it 
must be that the poles in e in their sum cancel. We 
have verified that this is indeed the case. Figs. nie,f) 
yield a pole in e which cancels against the second order 
contribution from Fig.^b) after the charge renormaliza- 
tion in (|2.6() is accounted for. The self-energy and vertex 
corrections in Figs, ^g-i) cancel against each other, as 
usual. And finally, Fig^j) represents an 'interference' 
between two first order contributions, the bulk renormal- 
ization of the vertex and the finite impurity contribution 
in Fig. ^a), and so requires no additional renormaliza- 
tions. Collecting these results, and using the fixed point 
value of / from 12.8|l . we obtain 



imp "imp 



727r2Q2 
7V2 



+ 0(e'), 



(2.27) 



although Aj'Jjjp and A^J^p are not expected to be equal at 
higher orders. 



C. Linear response to a uniform applied field 

We will be interested in the spatially dependent mag- 
netization response to an external static uniform mag- 
netic field Hu — H. Let us start with the calculation of 
the impurity susceptibility Ximp, which characterizes the 
total excess magnetization arising due to the impurity in 
response to the applied field. 



1. Impurity susceptibility 
The susceptibility can be defined as 



Xi 



imp 



7V2 



a 



(2.28) 



where the prefactor stems from averaging over the "spin" 
space of SU(N) generators. Using the action (|2.1(l . modi- 
fied according to H2.14|l , one can see that the presence of a 
uniform magnetic field leads to three new vertices shown 
in Fig.|2a-c). The impurity-related diagrams must con- 
tain at least one impurity "source" vertex shown in Fig. 

lae). 

It is easy to show that the leading order correction to 
the free energy is given by the diagrams of Fig. EJa-g). 
Notice that the leading contribution is of the second order 
in the impurity charge Q, since all first-order diagrams 
happen to be odd in Matsubara frequency and thus are 
identically zero. The propagator enters those dia- 
grams only as D^o{k,ULi = 0) = Sf^o/k"^, and the latter 
expression is valid for any choice of gauge which ensures 
that our final results will be gauge- independent. Sum- 
ming up the contributions of the diagrams (a-g) of Fig. 



13 one obtains 



d'^k dfik' 1 



2 



1 



16^2 16^4 



K+x2)2(c.2+y2) (^2 +^2)3(^2 +y2) 

- 2tr2 2. 1 (2-29) 

where = |fc|, and the sum in g{x, y) runs over bosonic 
Matsubara frequencies w„ = 27:111 . The overall factor of 
S comes from a trace over the components of the field z 
and is defined as 



c 



7V2 - 1' 



(2.30) 



where C is the eigenvalue of the Casimir operator of 
SU(N) which depends on the specific representation of 
the group. In the SU(2) case, this factor takes the famil- 
iar form S = S{S + l)/3. It is worth noting that all the 
diagrams are symmetric with respect to the exchange of 
k and k + k' , so only the symmetric part of g{x, y) yields 
a finite contribution to the susceptibility. 

Obviously, the expression (|2.29|) contains a singularity 
at k,k' ^ 0. At T — this singularity is formally unim- 
portant, since at T ^ the function g(x,y) becomes 
fully antisymmetric: 



\im^g{x,y) 



x-y 



2xy{x + yY 



(2.31) 



so at T = the impurity susceptibility is identically zero. 
Alternatively, one may notice that the expression under 
the sum sign in (|2.29|) is a full derivative in which 
implies that Ximp vanishes at T = when the sum in a;„ 
is replaced by an integral. 

At finite T, however, the symmetric part 'g{x, y) — 
1/2(5(2:, y) + g{y, x)) is finite: 



gix,y) 



a;(3y2 -f a;2)cosh(y/2T) 



+ {x ^ y), 



16xy{x^ - y^)T^ sinh^{y/2T) 

(2.32) 

and the infrared singularity becomes dangerous (the in- 
tegrand in behaves as k~^{k')-^ at k,k' 0). 
Such finite-T infrared divergencies in massless theories 
are well-known^ and can be cured by incorporating the 
thermal "screening mass" into free propagators. The 
thermal masses nip and mf, (for the "photon" and 
for the boson z, respectively) can be easily found to the 



(a) 



(b) 




FIG. 4: Self-energy diagrams determining the thermal screen- 
ing masses (12.3311 for the "photon" (a) and for z-boson (b). 



first order in the coupling constants Bq and ug by comput- 
ing the temperature-dependent contribution to the self- 
energy diagrams shown in Fig. ^ 



6 



-T2, 



|_(A^ + l)uo 



12 



^}r2. (2.33) 



At the critical point Cg and Uq are both proportional to 
e = 4 — D, so that we can write 



(2.34) 



where Ap.b can be found from the beta- functions (|2.8|1 : 
one obtains Ap = 47r2 , and the expression for Af, is rather 
complicated'*2 but simplifies in the large- limit: 



An 



47r2, 



Ab -> 7r2/6, N 



(2.35) 



Taking into account the thermal masses amounts to 



replacing {k' 



{k'Y 



li and 



9 9 
ml, y^ 



y + ra^ in Eq. (|2.29|) . It is easy to see that the main 
contribution to the integral in H2.29|l comes from small 
momenta k,Q < p <g; T, so one can replace ]){x, y) by 



lim o(a;,2/) 



1/2 



1/2 



{x^ 



m2)2 



+ 



(a;2 -I- ml){y^ + ml) ' 



(2.36) 



This means that the main contribution to the suscepti- 
bility comes from the zero Matsubara frequency {ujn = 0) 
part of the expression H2.29|l and thus is determined by 
the first two diagrams in Fig. O Evaluating the momen- 
tum integrals in H2.29(l in c? = 3, and substituting Cq by 
its fixed-point value 2^^ ^ one readily obtains the Curie- 
like law (|1.3|) for the impurity susceptibility, with the 
anomalous Curie constant C given by 



C = 



9eQ25'A2 
27r2Ar2 



1 



■ \ + 2-\/ApAfc 4-y/ ApAfc 



}. (2.37) 



It is worthwhile to note that the resulting expression for 
Ximp is proportional to e and not to as is apparent 
from (|2.29() . because of the partial cancellation caused 
by e-dependence of the thermal masses (|2.34|l . 



2. The Knight shift 

The magnetization response to a uniform magnetic 
field is space-dependent due to the presence of an impu- 
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rity; experimentally, this can be detected in NMR exper- 
iments by measuring the Knight shift. It is not dilRcult 
to generalize the calculation of the previous subsection 
to the space-dependent case; the "uniform" (slowly vary- 
ing on the scale of a lattice constant) component of the 
Knight shift Ku{x) is given by 



Ku{x) = / d''x'xu{x,x'), 



where Xu{x, x') is a generalization of H2.28|l . 



Xu{x,x') = 



T 



N2 - 



<52 InZ 



SHuaix)SHua{x')' 



(2.38) 



(2.39) 



This uniform Knight shift obeys the scaling form in 
Eq. H1.4|l . A direct application of Eq. H2.22|l to this scal- 
ing form implies, as in Ref.l^that in the quantum-critical 
region 



Ku{x ^ 0) 



1 



2-l-»;'/2 |^|<i-7,'/2 ■ 



(2.40) 



One can see that Xuix,x') is determined by the same 
set of diagrams of Fig. Ela-g) but with momenta flow- 
ing into the external vertices. To obtain the Knight shift 
from Eq. (|2.38|) . it is useful to define the Fourier trans- 
form Xu{q) by 



Kuix) 



(27r)" 



Xuiq)e 



iqx 



(2.41) 



and then the diagrams have momentum q flowing out of 
one external vertex. Obviously, Ximp ~ 

(g — 0) and 

the normalization (|1.5() is satisfied. Similarly to the case 
of susceptibility considered above, one can show that the 
Knight shift is identically zero at T = 0. At finite tem- 
perature the main contribution to Xuifl) is given by the 
diagrams of Fig. |2Ia,b) with the Matsubara frequency in 
each propagator replaced by the corresponding thermal 
mass (|2.34(l . yielding 



Xu{q) - 2Q\ef^^ST[l{q,mp)I{q,m^) (2.42) 

I{k,mb) I 



{2TTy (fc2 + m2)[(fc - q)2 + 

where we have denoted 

d'^k 1 



I{q,m) 



{2itY {k^ + m2)[(fc + g)2 + m^] ' 
We will only need the expression for I{q, m) in rf = 3: 

/fc3 (<?,«) = arctan M . (2.43) 

47r|(7| 2m 

Asymptotic expressions for Xn{q) can be easily obtained. 
For small wave vectors q <^ 2mi, one recovers the results 



for the impurity susceptibility from the previous subsec- 
tion, Xu{q *C 2m{,) ~ Ximp, while for large wave vectors 
q ^ rrip one has 



Xuiq) ^ ST 



12TTQe\ In^q/iTip 



N 



q- 



q > mp. (2.44) 



At large distances r ^ ^ the Knight shift decays ex- 
ponentially as e"™'"'' due to the finite photon mass nip. 
At small distances the Knight shift has the following be- 
havior: 



K^{r)^367rST[^] -\n(— 
^ N / r \ 7n„r 



riipr <€. 1. 



(2.45) 

Notice that nip oc eT^, so at low temperature the condi- 
tion nipr <^ 1 will be valid for a wide range of r. Apart 
from the log correction, (|2.45|) is consistent with the gen- 
eral scaling forms H1.4|) . (|2.22|) . and H2.40|l : the logarithm 
is likely an artifact of the integer-valued exponents that 
appear at low orders in the computation. 

In a similar manner, one can calculate the staggered 
Knight shift 



Ks{x) 



^2 InZ 



5HuaSHsa{x)' 

a ^ ' 

This obeys a scaling form analogous to Eq. (|1.4|l 



(2.46) 



where the behavior as x 



rT^(D-2+r,)/2 

Ks{x) $.(Tx) (2.47) 

analogous to Eq. (|2.40|l is 
1 1 



(2.48) 



The staggered Knight shift is computed by the Fourier 
transform of the momentum-dependent staggered suscep- 
tibility Xs{q) which is in the leading order determined 
by the diagrams of Fig. |2th,i), with an uncompensated 
momentum q being "injected" in the staggered magnetic 
field vertex of Fig. Ofd). Notice that for the staggered 
susceptibility the first nonvanishing contribution appears 
already in the first order in the impurity charge Q. As 
in the uniform case, the staggered susceptibility can be 
proven to be identically zero at T = 0, and at finite 
temperature the result is again determined by the zero 
Matsubara frequency term and thus by the diagram Fig. 
|3fi) , with the zero frequency replaced in each propagator 
by the respective thermal mass. This yields 



Xs{q) 



ien^QST I{q,mb) 
N q^ + rni 



(2.49) 



Using (|2.43|l . one can easily deduce the asymptotic be- 
havior of the staggered Knight shift at small distances 



r <^ m 



p 



3eQST 
Ks{r) = ^ In 



mpT <^ 1. 



(2.50) 
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Again, apart from the logarithmic correction, this result 
is in accord with the scaling form suggested by H2.21|l and 
(E3H1)- It is worthwhile to note that in CP^"^ model for 
small e = 4 — D and large N the staggered magnetiza- 
tion response Kgix) to an external magnetic field is much 
stronger than the uniform response Ku{x). 



D. Neel ordered phase 

In the Neel phase, one of the components of the z field 
gets condensed. Assuming (zq) oc 5ai, we can parame- 
terize 



^J\s\/u + cr j e'^ for a = 1 
ipaB^'^ for a > 1 



(2.51) 



where cr is a real field, and ip^, ct = 2 . . . N, are — 1 
complex fields. The overall phase C can be gauged away, 
and we can set C = 1- Then becomes a massive gauge 
boson, with the "mass" \/2e^|s|7u acquired by the Higgs 
mechanism. Inserting (|2.51|) in Sb -f <Simp we can obtain 
a straightforward perturbative expansion for all physical 
properties in powers of uq and Bq. 

First we describe the behavior of the Neel order near 
the impurity. We expect this to obey the scaling form 



(2.52) 



where (3 = {D — 2 + rj)v /2. For large x, we expect that 
approaches a constant which characterizes the bulk order, 
while for small x near the impurity we may deduce from 
(IT^ that 



1 



y 



(D-2+r^-q')/2 



(2.53) 



Explicitly, computing {<j)a) by the perturbative expansion 
defined above, we find that there are no impurity depen- 
dent corrections at first order beyond tree level. So to 
this order, F^{y) is independent of ?/, but y dependent 
terms do appear at higher orders in e. The spatially in- 
dependent response at this order is also consistent with 
the exponent relation H2.19|l and (|2.53|) . 

Turning next to the uniform magnetization, we now 
expect from the scaling relations in Section III BI that 



{Ma{x)) = \s\''-'F„,{x/\sn 
where now the scaling function obeys 

\imFra{y) 



d-ri'/2 



(2.54) 



(2.55) 



It is easy to compute Ma to one-loop order, and we find 



(Maix)) 



QTfi2e2|s| f d'^q 



piqx 



{2nY q' + 2\s\ef,/u 



(2.56) 



This obeys the scaling form in (|2.54|) with a magneti- 
zation that decays exponentially on a scale set by the 



"photon" mass. For small x, note that to leading order 
in e, {Ma{x)) ^ and this is consistent with (|2.55|) 
and the value of 77'. 

An interesting property of (|2.5t)|) is that the total mag- 
netization near the impurity is given by 



(2.57) 



As in Ref. ,9, this relation is expected to be exact, and 
the total magnetization near the impurity is therefore 
quantized. The quantization is a consequence of the con- 
servation of total spin. Here it may be viewed as an inter- 
esting consequence of Meissner flux expulsion induced by 
the condensation of the Za ■ In the presence of an applied 
field Ha, the fact that the Za condensate is polarized 
along the a = 1 direction, the Meissner effect implies 
from (|2.14|l that iAr + HaT^i will fiuctuate around zero. 
So setting Ar = iHaTf^, we obtain from 1)1.2(1 a term 
in the effective action which is linear in Ha, and whose 
coefficient is the total magnetization in H2.57|) . 



E. Spin gap phase 

Finally, we turn to the spin gap phase which appears 
for sufficiently large s. We will not include 'danger- 
ously irrelevant' monopole effects hereipSiSSiSi which are 
believed to lead to VBS order and confinement at the 
longest length scales. We will restrict our attention to 
length scales of order 1/A, where A is the spin gap. We 
will argue that in d = 2 a spinon is necessarily bound to 
the impurity, and so the onset of confinement at longer 
scales does not have a significant infiuence on the behav- 
ior of the impurity. 

It is straightforward to extend the e expansion to the 
spin gap phase. The spinous experience the Coulomb 
potential of the impurity ^ This potential may 

form bound state, but because of the bulk spin gap, these 
states are above the ground state. So the spin gap is pre- 
served in the presence of the impurity and all magnetic 
response functions vanish exponentially at low tempera- 
ture. 

We now argue that these conclusions of the e expansion 
are strongly modified in c? = 2. The key point is that the 
Coulomb potential ^ ln(r) in d = 2, and this increases 
without bound as r — )■ cxd. Consequently the self energy 
of a naked impurity with nonzero charge diverges loga- 
rithmically with system size. Therefore, it always pays 
to create a spinon above the spin gap, and bind it to the 
impurity, and neutralize its charge. The scale at which 
this occurs can be estimated using the 1/A expansion4^ 
In this expansion, the Coulomb interaction is screened 
by the vacuum polarization of spin-anti-spinon pairs to 
a value of order (QA/A) ln(i?A). The spinous obey a 
Schro dinger equation in this potential with a mass A, 
and from this we obtain an estimate for the size of the 
bound state of order R ~ {N/Qy/'^/A. At length scales 
larger than this R, the impurity will just appear as an 
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isolated spin moment, and will therefore contribute the 
ordinary Curie susceptibility = 1/(4T) (for a Q = ±1 
impurity) . 

Thus the impurity simply behaves as a free local mo- 
ment in the spin gap phase. Note that we did not have to 
appeal to confinement physics to reach this conclusion; 
confinement, and VBS order, appears at the paramet- 
rically larger scale^- of order a(a/A)"^, where a is the 
short-distance cutoff. 



III. STAGGERED FLUX SPIN LIQUID 

This section will examine the response of a second al- 
gebraic spin liquid to impurities: the staggered fiux spin 
liquid with fermionic spinon excitations. ^^ '^^•■^^■^^ We will 
see that the scaling structure is rather similar to that 
of the critical CP^~^ model examined in Section^ al- 
though there will be some crucial differences in some 
physical properties. 

The low energy excitations of the spin liquid are de- 
scribed by Nf flavors of 2-component Dirac fcrmions, ^1/, 
coupled to the U(l) gauge field with the action 



(3.1) 



As before, y = (r, x) is the spacetime coordinate, and fi 
extends over the D = 3 spacetime indices. The Dirac 
matrices 7^ — (t^,t^, — t^), where r'^ are the Pauli ma- 
trices, and the field ^ = V"^ — we follow the notation 
of Hermele et al^ The number of flavors is Nf = 4 for 
the staggered fiux state, and our results are also extended 
to the so-called 7r-flux state, which has Nf = 8. Our 
analysis will be carried out, as in the previous works, in 
a 1 /Nf expansion. We have not included a bare Maxwell 
term for the gauge field in iS;, because it turns out to be 
irrelevant at all orders in the 1 /Nf expansion. 



A. Bulk theory 

The structure of th e bulk theory has been described 
in some detail in Refs. I35l37l and we will not repeat the 
results here. In the large Nf limit, the propagator of the 
gauge field in the Lorentz gauge is 



DM 



p2 



16 



(3.2) 



This propagator arises from the vacuum polarization of 
the fermions. Notice that it is suppressed by a power of 
l/Nf, so the 1/Nf expansion can be setup as a pertur- 
bation theory in the fermion-gauge field interaction. 

A large number of order parameters can be constructed 
out of fermion bilinears, and a detailed catalog has been 
presentedj^ Among these are the SU(A^/) fiavor currents 



ja 



(3.3) 



certain components of which are the magnetization den- 
sity and current. Conservation of this current implies the 
scaling dimension 



dim[J^] = 2, 

which is the analog of the relation (|2.1t)|) in d ~ 
will restrict all discussion in this section to d = 2) 
Also considered were the following quantities: 



N" 



(3.4) 

2 (we 

(3.5) 



which relate to additional order parameters including the 
Ncel order (note that AI does not correspond to the phys- 
ical magnetization). Their scaling dimensions have been 
computed in the 1/A^/ expansion: 



dim[A^"] = 2 



64 



37r^A^/ 



0{l/N'f) 



(3.6) 



B. Relationship to earher work 

Before describing the results of our analysis of the im- 
purity in the U(l) sF phase, it is useful to describe the 
earlier analyses in Refs. 24 25-26 27, .2&. They ignored 
the Af^ fiuctuations, but instead considered a theory of 
fermionic spinous fia on the sites, i, of the square lattice; 
a is a spin index. These spinous obey J2aifiafia) = 1 
on every lattice site except at the impurity i — 0. Here, 
they inserted a vacancy (representing e.g. a Zn impurity) 
by including a potential term in the Hamiltonian 



imp 



(3.7) 



and taking the limit ^ cxd to prohibit any spinous from 
residing on the vacancy. The key physical ingredient in 
these analyses is the difference in the spinon occupation 
number between the impurity and the bulk: 



^ ] { fia fia) 



1 



JiO- 



(3.8) 



In our approach the analog of the V ^ 00 limit is ob- 
tained by the functional integral over Ar . In the bulk the- 
ory, the continuum field 'J is defined so that ('5 ^'5) = 
in the absence of the impurity, and the continuum analog 
of Eq. 1231) is 



(^^^(r)) = -QS'^ir) 



(3.9) 



where the rhs has a Dirac delta function. The constraint 
in Eq. I|3.9|) is imposed not by adding a potential energy, 
but by the functional integral over Ar, which appears 
in the Lagrangian density of Sb + Simp (with Sb given in 
Eq. ITTJ) a,siAr{^''^+Q6'^{r)). Our treatment of gauge 
fluctuations ensures that the constraint on the spinon 
occupations is imposed not just on the average, but dy- 
namically on all states and on all sites. The equality in 
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FIG. 5: Feynman diagrams to order 1/Nf for the Ihs of 
Eq. 13.911 . The fuh Une is the Dirac fermion propagator, the 
wave line is the photon propagator, the filled square is the im- 
purity source term in <Simp, and the filled circle is ^'^\[' vertex. 
As noted in the text, the diagram in (b) has an odd number of 
photon vertices, and so vanishes by Furry's theorem.** Here, 
and henceforth, we do not show a number of other diagrams 
which vanish because of Furry's theorem. 



Eq. H3.9|) holds to all orders in the 1/A^/ expansion: this 
follows from the 'equation of motion' for A^- 



^ 0. 



(3.10) 



It is instructive to also test Eq. (|3.9|) by explicitly eval- 
uating the Ihs of Eq. (|3.9|) in the \/Nf expansion. The 
corresponding Feynman diagrams are shown in Fig. 
At leading order, we have the diagram FigEfa), which is 
easily evaluated to yield the rhs of Eq. (|3.9I) . At order 
1/Nf, the diagrams shown in Fig.|3Ib-g) contribute. Of 
these, Fig.lSJb) vanishes by Furry's theorem^i Of the re- 
maining, it is easy to show that they cancel in pairs: this 
requires only the knowledge that the photon propagator 
is the inverse of the fermion vacuum polarization bubble. 
Thus Figs. Etc) and (d), (e) and (f), and (g) and (h), 
all cancel against each other, and Eq. (|3.9II is thus es- 
tablished to this order. It is not difficult to extend these 
arguments to all orders in l/Nf. Note that it is possible 
to satisfy Eq. (|3.9|) in a perturbative treatment of gauge 
fluctuations, without the need to appeal to bound states: 
the delta function at r = arises from a superposition 
of the contribution of many extended states, rather than 
from a bound state claimed earlier 

It is clear that our approach yields a systematic 
and controlled treatment of the spinon deficit at the 
impurity, in contrast to the earlier ad-hoc mean-field 
approaches i2ii2S*2SiSL2& Our analysis implies that the 
Curie constant C is a non-trivial number, with contri- 
butions at all orders in 1/Nf, and is not given simply by 



the Curie constant of a single spin. 

The proper analysis of a potential term like that in 
Eq. (|3.7|l requires the scaling analysis of perturbations 
to the conformal field theory defined by Sb + Simp- The 
action for such a perturbation takes the form 



V / dT*^*(2; = 0,r) 



(3.11) 



A simple analysis of scaling dimensions shows that 



dim[y] = -l + C'(l/iV/). 



(3.12) 



So, potential scattering is an irrelevant perturbation at 
all orders in the 1/Nf expansion. 

A further distinction between our results and the ear- 
lier ■vvorli2li2SiSS42L2& appears in the response to a uniform 
magnetic field. This has significant experimental conse- 
quences, and will be discussed in Section IlII Dl 



C. Impurity exponents 

The analysis is analogous to that carried out in Sec- 
tion lllBl for the CP^^^ model. As there, and for similar 
reasons, we neglect the possible influence of monopoles. 

Here, we will determine the impurity renormalization 
to the order 1/Nf. To this order, we find that all the 
operators introduced above, the and N'^ acquire a 
common correction, Ajmp, to their bulk scaling dimen- 
sion. This correction is given by the single diagram in 
Fig n^c) , which (along with its symmetry related part- 
ner) evaluates to 



512Q2 



d'^k {—iLtj + q-f){-iuj+{k + q)-f) 



(27r)2 (27r)2 



kq{q'^ + cj2)((q + k)^ + uj^) 



(3.13) 



where r = (r^, r^) and A is an ultraviolet cutoff; here and 
henceforth we are using the fermion lines to represent the 
propagator (iu) + t ■ q)^^ from ^' to ^'^ (rather then the 
Dirac propagator from 5* to ^P). There are a number of 
other diagrams, like those shown in Fig. ^ which could 
contribute to the vertex renormalization; however they 
do not contribute either for reasons as discussed for the 
(^pW-i niodel, or because of Furry's theorem. From this 
we obtain the impurity correction to the scaling dimen- 
sion 



A, 



128Q2 
Nj^ 



0{1/Nf) 



(3.14) 



We expect that the higher order corrections will not be 



the same for the J° and N°-. 
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FIG. 6: Feynman diagrams for the impurity susceptibility 
and Knight shifts of the staggered flux spin liquid. The grey 
circle vertices depend upon the correlator being evaluated: 
they equal (i) T'^'yr for the ferromagnetic spin density J", 
[ii) T"' for the order parameter A''", and (iii) unity for M. 



D. Linear response to a uniform applied field 



Again, we proceed in analogy to the analysis in Sec- 
tion We apply a uniform magnetic field Ha, as- 
sociated with the SU(A^/) generator T", which couples 
linearly to the conserved total spin density J". In prin- 
ciple, in the presence of the impurity, the linear response 
to this applied field can induce time-independent, space- 
dependent average values of not only the spin density, 
J!^, but also of the other order parameters M, and A^°. 
This would be analogous to the non-zero averages of the 
uniform and staggered magnetizations induced by an ap- 
plied magnetic field on the CP^~^ model (which led to 
the uniform and staggered Knight shifts). However, here 
we will find a crucial difference. In the scaling limit of the 
algebraic spin liquid represented by Eq. (|3.1(l , Ha induces 
only a non-zero J", and average values of all the N°- and 
M are zero. Thus there is only a uniform Knight shift, 
Ku{x), and all the 'staggered' Knight shifts, Ks{x), as- 
sociated with the many competing orders are zero. The 
staggered Knight shift can appear only if some correc- 
tions to scaling are included, associated with irrelevant 
operators which reduce the symmetry of the conformal 
theory to that of the lattice model, and so can be ex- 
pected to be weaker than the uniform Knight shift. 

The vanishing of {M{x)) and {N°'{x)) in the presence 
of Ha can be established by a careful consideration of 
the symmetries of the Dirac fermion theory. First (Af ) = 
0, to linear order in Ha, simply by SU(A^/) symmetry. 
Establishing the value of {N"") requires more complicated 
considerations. Let us consider the leading contribution 
to (iV°), to linear order in Ha, in the 1/iVj expansion, 
represented by the graph in FigEJa). The value of the 



fermion loop is proportional to 

/ ^Tr[(*c^„ + f.fc)-V3(ic^„+T.fc)-i 
x{iujn + f ■ {k + q))-^] (3.15) 

Evaluating the trace over the Dirac matrices, we obtain 
an identical zero. This can be understood as a conse- 
quence of time-reversal invariance. Both Ha and iV° are 
odd under the time reversal;^ as is the charge of the 
impurity Q. So an expansion of {N°')/Ha can only in- 
volve even powers of Q. Proceeding to order ^/N'j, we 
obtain diagrams like those shown in Fig. |B{b-d), which 
have a pre- factor of , and so are potentially non-zero. 
However, these diagrams have a fermion loop with an 
odd number of 7'^ vertices, and so vanish by Furry's 
theorem44 By a combination of Furry's theorem and 
time-reversal invariance we can now easily see that all 
terms vanish and so (A^°) = 0. Because of the ma- 
trices in the definitions of iV" and J°, there must be a 
single fermion loop which connects the external vertices. 
By Furry's theorem, this loop must have an odd number 
of photon vertices. All other fermion loops can only have 
an even number of photon vertices. Consequently, there 
must be an odd number of photon vertices remaining to 
connect to the external impurity source term. However, 
by time-reversal, there must be an even number of impu- 
rity terms; hence the result. 

Our conclusion that (iV") = is starkly different from 
the mean-field prediction2li2Si2Si2LS£ of an induced mo- 
ment which oscillated strongly between the two square 
sublattices. Such oscillations can only appear upon in- 
cluding irrelevant operators. 

It remains only to compute the uniform Knight shift 
Ku{x), or equivalently, its Fourier transform x«('?) de- 
fined in Eq. H2.41|l . The scaling analysis of Section llll CI 
implies that this Knight shift obeys the scaling form in 
Eq. H1.4|l . with the x ^ behavior given by 

Ka{x ^ 0) - Td-i+iy-,^ 1 (3^^g) 

\x\ ^'"P 

Here we expect that the Ai^p exponent is that associated 
with J", and not (unlike the situation with the CP^^^ 
model) that associated with the N°' , because there is no 
'mixing' between the staggered and uniform magnetiza- 
tions near the impurity for the sF spin liquid. 

At first order in 1/Nf, Xu{q) is given by the diagram in 
Figinia), which vanishes by Furry's theorem. The leading 
non- vanishing contribution is at order 1/N'j, and is given 
by the diagrams in Fig|HIb-d) (a number of order ^/N'j 
diagrams which vanish because of Furry's theorem are 
not shown). The sum of these diagrams can we written 
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in the following compact form 



d?k (Pp d 



-Tr 



47r2 47r2 i9(iw„) 

xDrrip)Drri\q- P\)- 



(3.17) 



Here Dtt is the r component of the photon propagator, 
representing the external lines connecting to the impu- 
rities, and S is the constant proportional to the SU(N) 
Casimir eigenvalue as defined in H2.30|l . Notice that the 
expression is a total frequency derivative; this immedi- 
ately implies that Xu{q) = at T 0, when the fre- 
quency summation can be converted to an integration. 
This vanishing is a consequence of the conservation of 
total spin. 

A non-zero Xuil) is obtained at T > 0, and we now 
compute this. First we need the photon propagator at 
T > 0. This is given by a single fermion loop and at 
this order we only need the t, t component at a spatial 
momentum q: 



iojn + f ■ k 



X (^iuj„ + T-{k + q}^ 



(3.18) 



We first combine the denominators in Eq. H3.18|l using 
the Feynman parameter u, perform the frequency sum- 
mation, and finally integrate over k. This yields 



NfT 



du In 



2 cosh 



(3.19) 

Eq. H3.19I) interpolates between Njq/IQ for T ^ q which 
agrees with the T = result in Eq. to (r/7r)ln2 

for g < T. 

Inserting Eq. (|3.19l) into Eq. I|3.17|l . we conclude that 
Xu{,q) obeys the scaling form 



Xu{q) = ^<^u[q/T) 



(3.20) 



where the scaling function $„ is the Fourier transform 
of the scaling function in Eq. H1.4|l . Further, C, the 
anomalous Curie constant appearing in Eq. (|1.3(l , equals 
$„(0). To determine the scaling function $„, we need to 
evaluate Eq. H3.17|l . We combined the Green's functions 
using two Feynman parameters u, v, evaluated the inte- 
gral over k, differentiated with respect to frequency, and 




FIG. 7: Numerically calculated scaling function $ii(y) de- 
fined by 13.1913.231 . The inset shows large-j/ behavior to- 
gether with the least square fit (dashed line) to a power law 
for y > 40: the fit yields 5„(?/) ~ C /y" with C = 24.7(1) 
and Q = 2.070(5). 



finally evaluated the summation over ojn- This yields 



Xuiq) = -^-^ f du f dv 



TT 



"'0 
2 „,„2 „,„2 



/'(A) 



47r2 1 

2A2 - uq'^ - vp^ - (l - 2(m + v))p- q 



4A 



>Drr{p)Drr{\q-'P\) (3.21) 



where 



q^u{l -u)+ p^v{l -v)- 2uvp- q (3.22) 



and / is the Fermi function 

f{e) = l/{e^'^ + l). 



(3.23) 



The remaining integrals have to be evaluated numerically. 
From such an evaluation aX q — Q we found 



1.0460(5)S'g2 



0{llN}) 



(3.24) 



for the universal Curie constant appearing in Eq. p.3|l . 
The calculated shape of the scaling function $u(2/) is 
shown in Fig. [T] From the numerical analysis of its be- 
havior for large arguments, we deduce that <&«(?/) has a 
power-law decay (y) fx 1/ at y 3> 1 , where the ex- 
ponent a K, 2 within the accuracy we were able to achieve 
when calculating the four-dimensional integral in l|3.21|l . 
Therefore, we obtain that at small distances the uniform 
Knight shift Ku{x) tends to a (lattice cutoff-dependent) 
constant as x — > 0, which is consistent, to the leading 
order in 1/Nf, with Eq. H3.16|l and the fact that Ajmp 
vanishes at this order (see Eq. (|3.14|) '). 
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IV. OTHER SPIN LIQUIDS 

This section will discuss the extension of our results to 
some other examples of spin liquids. The cases consid- 
ered below do not share the high degree of universality 
found in the critical CP^~^ and sF spin liquids above: 
some aspects depends upon microscopic details. We will 
consider the spinon Fermi surface in Section flV Al and Z2 
spin liquids in Section ITVBI 



A. Spinon Fermi surface 

This spin liquid differs from the sF spin liquid by hav- 
ing a Fermi surface of gapless spinon excitations, with 
a non-zero density of states. The bulk theory is not a 
CFT, and the bulk action for the fermionic spinous is 



Si 



dr 



{2n) 



+e{k-A) Uik,T) (4.1) 



Here a is a SU(7V) spinor index, e(A:) is the spinon dis- 
persion, and the Fermi surface defined by e{k) = has 
a finite density of zero energy spinon excitations. The 
properties of this spin liquid can be computed by a com- 
binations of previous methods^^*^ and the methods in- 
troduced here. 

In the absence of the impurity, the spin susceptibility 
is finite as T ^ 0, equal to the density of states at the 
Fermi level. The impurity introduces a 'Coulomb' po- 
tential for the spinous. In the small momentum limit, 
k ^ 0, this interaction is screened to a constant value by 
the Fermi surface excitations. Evaluating the response 
to an applied magnetic field using a diagram like that 
in Fig. EJa), we conclude that Ximp(?' 0) is finite and 
proportional to the energy derivative of the density of 
states at the Fermi level. 

Upon considering the spatial dependence of the Knight 
shift, again using the diagram in Fig. Efa), we find 
that the dominant response oscillates at the '2kp^ 
wavevectorn^ The transverse gauge interactions strongly 
enhance the impurity-spinon vertex, and this leads to a 
singular response at 2kp d ecaying with a power-law of x 
determined in Refs. I45l4(l Again, this response is finite 
as T 0, although the singular portion will have strong 
corrections at T > 0. 



B. Z2 spin liquid 

We will describe the Z2 spin liquid accessed from the 
(^pA^-i jy^odel of Section |n] by condensing a charge-2 
Higgs scalar, ^Mdl We add a term like \{d^ - 2zA^)$p 
to the full action of the CP^~^ model, and assume we are 
deep in the phase with ($) ^ 0. This condensation gaps 
out the Af^ excitations, and so it is useful to completely 



integrate out the A^. The result is an effective action 
for bosonic spinons Za which describes a bulk transition 
from a ordered state with non-collinear magnetic order to 
a Z2 spin liquid. This is proposed as a possible model for 
the triangular lattice antiferromagnetii2i^2iffi The effec- 
tive action obtained by this method has the bulk action 

S,^J d^y{\d,z^\^ + sK\' + '-f (|z„|2)' 

+Vi K^^ZqI^ -I- V2 {z^d^Zaf + c.c.j, (4.2) 

and an impurity term which descends from its U(l) 
charge 



Simp — I dr 



-|Aiz* (a; = 0, T)drZa{x — 0, r) 

+A2|za(x = 0,r)|2}. (4.3) 

Here the coefficient Ai is proportional to the impurity 
charge Q, but the value of Ai is not quantized. An initial 
derivation does not immediately yield the A2 term, but 
it is clearly allowed by the symmetries and is evidently 
more relevant than the Ai term: it represents a poten- 
tial scattering term for the bulk spinons. The critical 
point of the bulk theory alone is a CFT, and it is ex- 
pected that the vi^2 are irrelevant at this critical point, 
leading to global 0(2A^) symmetry'^^'^i. However, as we 
will see below, the perturbations in iSimp are not exactly 
marginal. 

It is useful to compare this theory with that of the 
recent work of Florens et al}^ Unlike us, they include an 
explicit spin degree of freedom at the impurity. In our 
approach, such a situation would be appropriate when 
the couplings Ai^2 are large enough to bind a spinon in 
the ground state at the impurity. In a Z2 spin liquid, 
there are no long range gauge forces (the Z2 gauge forces 
are short range), unlike the situation in Section [ll El and 
so it is not required that the ground state have non-zero 
spin for a vacancy impurity, e.g., Zn on a Cu site. For 
a magnetic impurity, such as Ni on a Cu site, a finite 
spin ground state may be more likely, but is also not 
required when the bulk system is in a gapped Z2 spin 
liquid. We will restrict our discussion here to the small 
Ai_2 region of (|4. 21) + 14. 31) . and this is equivalent to being 
in the Kondo-screcned phases of Florens et alm^. 

Let us now describe the phases of iSf,+5imp as the value 
of the tuning parameter is scanned. 

For s 3> 0, we have the gapped spin liquid phase in 
the bulk. The bulk spinons will experience a potential 
from the Ai^2 terms, and in the attractive case in c? = 2, 
such a potential always has a bound state. However, the 
binding energy can be less than the bulk spin gap, and 
hence the ground state remains a singlet and the spin- 
gap is preserved. All magnetic response functions will 
therefore be exponentially small at low T. 

At the bulk critical point, we need to examine the 
scaling dimensions of the perturbations in iSimp- In the 
e = 4 — D expansion, simple power-counting shows that 
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both Xi 2 are irrelevant. Consequently, we do not expect 
a Curie spin susceptibility, but a suppressed response 
determined by the scaling dimensions of Ai.2- In the 
D = 2 + e expansioujiSi we work with the constraint 
= 1, and hence the A2 term disappears. The Ai 
term is marginal at the tree level, and a more complete 
renormalization group analysis is needed to understand 
the flow of A; this, however, is beyond the scope of the 
present paper. 

Finally, for s ^ 0, in the magnetically ordered phase 
with coUinear order, we can carry out an analysis of the 
magnetic properties as in Section IIIDI We apply the 
ansatz in Eq. (|2.51|) to St + 5imp above, and compute 
all physical properties as a perturbation series in uq, vi^2 
and A1.2. The result shows that a spatial dependence of 
the non-coUinear order and the uniform (ferromagnetic) 
magnetization is induced near the impurity by Ai^2- The 
uniform magnetization requires the 'particle-hole' sym- 
metry breaking term Ai, and is non-zero even in the ab- 
sence of an applied magnetic field. However, unlike the 
situation in Section III Dl this ferromagnetic moment is 
not quantized, and takes a value dependent upon the 
bare values of the couplings: this is because spin rota- 
tion symmetry is completely broken in the non-coUinear 
ordered phasei^ 



V. CONCLUSIONS 

Let us summarize the basic physical characteristics of 
the impurity response of the various spin liquid states 
considered in this paper: 

(i) Neel-VBS transition: The impurity charac- 
teristics across this transition are qualitatively very 
similar to those found for impurities in confining 
antiferromagnetSfiiSiiSi although the underlying field the- 
ory is quite different, as are all the exponents and scal- 
ing functions. The spin-gap phase has the conventional 
Curie response of a local moment at the impurity; for the 
confining antiferromagnets this local moment is induced 
by confinement physics, while for the Neel-VBS transi- 
tion it is induced by the logarithmic 'Coulomb' interac- 
tion associated with the U(l) gauge field. The Neel-VBS 
quantum critical point has an anomalous Curie response, 
with a Curie constant C specified in Eq. (|2.37l) to lead- 
ing order in the e = 4 — D expansion. The response of 
the Neel phase is as found in the earlier case: a spatial 
varying staggered and uniform moment is found near the 
impurity, with the moment in the latter quantized as in 
Eq. H2.57|l . The Knight shift has both staggered and uni- 
form components, with its behavior in the vicinity of the 
impurity specified by the impurity exponents determined 
in Section FlI Bl At the critical point, the induced stag- 
gered moment is much stronger than the uniform magne- 
tization component and thus represents the dominating 
response to an external field. 

(ii) Staggered flux spin liquid: This spin liquid is gener- 
ically critical, and so has an anomalous Curie response; 



the Curie constant is given, to the leading order in the 
1/Nf expansion, by Eq. H3.24|l . The Knight shift now has 
only a uniform component, but no staggered component 
in the scaling limit. The absence of any staggered Knight 
shift is a defining characteristic of the sF spin liquid, and 
is a consequence of its large emergent symmetry^ There 
is a large number of competing order parameters, and the 
leading impurity action iSimp has no natural way of choos- 
ing among them, leading to a response which is restricted 
to the ferromagnetic moment alone. The behavior of the 
ferromagnetic Knight shift upon approaching the impu- 
rity is specified by the impurity exponent l|3.14|l found in 
Section IlII CI to the leading order in 1/A^/, the Fourier 
transform of the ferromagnetic Knight shift is computed 
in Sect. lmPl 

(Hi) Spinon Fermi surface: This case has a finite mag- 
netic response as T — > 0, unlike the divergent Curie sus- 
ceptibility of the sF phase. The response decays with a 
power law away from the impurity, with the dominant 
response being at the 2k p wavevector. 

(iv) Z2 spin liquid: The spin liquid itself has an expo- 
nentially suppressed impurity response, because the spin 
gap is preserved in the presence of the impurity. The 
proximate ordered non-coUinear state has a response sim- 
ilar to that of the Neel state, but with no quantization 
of the total magnetization. 

Apart from their applications to spin liquids, the re- 
sults in Section IIIII also have a direct application to 
the physics of charged impurities in two-dimensional 
graphene. It is well known that the low energy elec- 
tronic excitations in graphene are described by 2 species 
of Dirac fermions. There is no fluctuating gauge fleld 
as in Eq. (|3.1() . However in the presence of charged 
impurity, the three-dimensional Coulomb potential the 
electrons experience has the form C /r, where C is a con- 
stant. Interestingly, this is exactly the form of the static 
potential found in Section IIIII in the presence of an im- 
purity, where {At) ~ 1/r. Thus, as long as we ignore 
quantum-electrodynamic loop corrections, the results of 
Section Hill applv also to graphene; speciflcally, in Fig. [SJ 
(a) and (b) apply to graphene while (c)-(h) do not. One 
of our important results for this system was Eq. I|3.9|l : 
that the induced charge density due to the Coulomb po- 
tential is a delta function, with Q a non-trivial univer- 
sal function of C. This results therefore applies also to 
graphene. It disagrees with the earlier result of Ref. IH^ . 
which found a 1/r^ decay in the induced charge density. 
We believe their results suffers from a cavalier treatment 
of the ultraviolet cutoff, which does not preserve gauge 
invariance. 

We now discuss possibilities for future theoretical 
work. 

For the sF and CP^^^ spin liquids, an alternative 
l/Nf (or respectively l/iV) expansion could give a useful 
complementary picture. We have used an expansion here 
in which the Nf 00 was taken at flxed Q. However, it 
is also possible to study the limit Nf — > 00 limit at flxed 
Q/Nf. The saddle-point equations are straightforward 
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to obtain for this case, but much more difficult to solve: 
they have the advantage of including the spinon-impurity 
'Coulomb' interaction at all orders already at Nj — oo, 
and so allow for spinon states below, or at the edge of, 
the bulk spinon continuum. We will report on the results 
of such an analysis in future work. 

It would also be useful to have a more complete under- 
standing of the s pin on Fermi surface case, by extending 
the work of Refs. l45ll46l to determine the gauge-invariant 
response functions near a single impurity. 

Finally, we comment on the experimental implications 
of our work. 

For the cuprates, NMR experimentsSSiSS show a large 
Knight shift response at the Neel wavevector in the vicin- 
ity of the impurity, and strong temperature-dependent 
enhancement of such correlations. These features are 
consistent with the Neel-VBS transition considered in 
Sectional but also with the transition in dimerized an- 
tiferromagnets considered earlier>i^'^^ The absence of a 
staggered response for the sF case is potentially a serious 



deficiency of the model of Section ITTll It remains to be 
seen if corrections to scaling (such as those in Eq. (|3.11|l ) 
can remedy the situation. 

For the organic Mott insulator k-(ET)2Cu2(CN)3, the 
NMR Knight shift'^^ has an appreciable T independent 
component at low T. This can potentially be fit either 
by the spinon Fermi surface, or by a weakly magnetically 
ordered state. 



Acknowledgments 

We thank A. Castro Neto, M. Hermele, Y. B. Kim, 
T. Senthil, and M. Vojta for useful discussions. This 
research was supported by the National Science Founda- 
tion under grant DMR-0537077. AK was supported by 
the grant K02335/1-1 under the Heisenberg Program of 
Deutsche Forschungsgemeinschaft. 



^ A. C. Hewson, The Kondo problem to Heavy Fermions, 
Cambridge University Press, Cambridge (1993). 

^ P. H. Nozieres and A. Blandin, J. Phys. (Paris) 41, 193 
(1980). 

^ B. A. Jones, C. M. Varma, and J. W. Wilkins, Phys. Rev. 

Lett. 61, 125 (1988). 
* L. Fritz and M. Vojta, Phys. Rev. B 70, 214427 (2004); 

M. Vojta and L. Fritz, Phys. Rev. B 70, 094502 (2004). 
^ C. L. Kane and M. P. A. Fisher, Phys. Rev. Lett. 68, 1220 

(1992). 

" S. Eggert and I. Affleck, Phys. Rev. B 46, 10866 (1992); 

Phys. Rev. Lett. 75, 934 (1995). 
^ S. Sachdev, C. Buragohain, and M. Vojta, Science 286, 

2479 (1999). 

® A. M. Sengupta, Phys. Rev. B 61, 4041 (2000). 

® M. Vojta, C. Buragohain, and S. Sachdev, Phys. Rev. B 

61, 15152 (2000). 
^° S. Sachdev, M. Troyer, and M. Vojta, Phys. Rev. Lett. 86, 

2617 (2001). 
" S. Sachdev, Physica C 357, 78 (2001). 
^2 L. Zhu and Q. Si, Phys. Rev. B 66, 024426 (2002). 

C. Zarand and E. Demler, Phys. Rev. B 66, 024427 (2002). 
^'^ Q. Si, S. Rabello, K. Ingersent, and J. L. Smith, Nature 

413, 804 (2001). 

M. Troyer, Prog. Theor. Phys. Supplement 145, 326 
(2002). 

" S. Sachdev, and M. Vojta, Phys. Rev. B 68, 064419 (2003). 
K. H. Hoglund and A. W. Sandvik, Phys. Rev. B 70, 
024406 (2004). 

S. Florens, L. Fritz, and M. Vojta, Phys. Rev. Lett 96, 
036601 (2006). 

T. Senthil, A. Vishwanath, L. Balents, S. Sachdev, and 

M. P. A. Fisher, Science 303, 1490 (2004). 

T. Senthil, L. Balents, S. Sachdev, A. Vishwanath, and 

M. P. A. Fisher, Phys. Rev. B 70, 144407 (2004). 

T. Senthil, L. Balents, S. Sachdev, A. Vishwanath, and 

M. P. A. Fisher, J. Phys. Soc. Japan 74, Suppl. 1 (2005). 

N. Read and S. Sachdev, Phys. Rev. Lett. 66, 1773 (1991). 



X.-G. Wen, Phys. Rev. B 44, 2664 (1991). 
S. A. Krivenko and G. G. KhaliuUin, JETP Lett. 62, 723 
(1995); G. KhaliuUin, R. Kilian, S. Krivenko, and P. Fulde, 
Phys. Rev. B 56, 11882 (1997). 

N. Nagaosa and T.-K. Ng, Phys. Rev. B 51, 15588 (1995). 
N. Nagaosa and P. A. Lee, Phys. Rev. Lett. 79, 3755 
(1997). 

C. Pepin and P. A. Lee, Phys. Rev. Lett. 81, 2779 (1998); 
Phys. Rev. B 63, 054502 (2001). 

Z. Wang and P. A. Lee, Phys. Rev. Lett. 89, 217002 (2002). 
J. Bobroff, H. AUoul, W. A. MacFarlane, P. Mendels, 
N. Blanchard, G. Collin, and J.-F. Marucco, Phys. Rev. 
Lett. 86, 4116 (2001). 
^° S. Ouazi, J. Bobroff, H. AUoul, M. Le Tacon, N. Blanchard, 
G. Collin, M. H. Julien, M. Horvatic, and C. Berthier, 
Phys. Rev. Lett. 96, 127005 (2006). 

A. Yazdani, B. A. Jones, C. P. Lutz, M. F. Crommie, and 

D. M. Eigler, Science 275, 1767 (1997). 

^2 S. H. Pan, E. W. Hudson, K. M. Lang, H. Eisaki, S. Uchida, 
J. C. Davis, Nature (London) 403, 746 (2000); E W. Hud- 
son, K. M. Lang, V. Madhavan, S. H. Pan, H. Eisaki, 
S. Uchida, and J. C. Davis, Nature (London) 411, 920 
(2001). 

Y. Shimizu, K. Miyagawa, K. Kanoda, M. Maesato, and 

G. Saito, Phys. Rev. B 73, 140407 (2006). 

S. Sachdev in Quantum magnetism, U. SchoUwock, 

J. Richter, D. J. J. Farnell and R. A. Bishop eds. 

Lecture Notes in Physics, Springer, Berlin (2004), 

cond-mat/0401041 

W. Rantner and X.-G. Wen, Phys. Rev. Lett. 86, 3871 
(2001); Phys. Rev. B 66, 144501 (2002). 

^® M. Hermele, T. SenthU, M. P. A. Fisher, P. A. Lee, N. Na- 
gaosa, and X.-G. Wen, Phys. Rev. B 70, 214437 (2004). 
M. Hermele, T. SenthU, and M. P. A. Fisher, Phys. Rev. 
B 72, 104404 (2005). 

^® F. S. Nogueira and H. Kleinert, Phys. Rev. Lett. 95, 
176406 (2005). 

B. L Halperin, T. C. Lubensky, and S.-k. Ma, Phys. Rev. 



17 



Lett. 32, 292 (1974). 

A. V. Chubukov, S. Sachdev, and J. Ye, Phys. Rev. B 49, 
11919 (1994). 

see, e.g., J. O. Andersen and M. Strickland, e-print 

hep-ph/0404164 , and refereces therein. 

I. D. Lawrie, Nucl. Phys. B 200, 1 (1982). 

G. Murthy and S. Sachdev, NucL Phys. B 344, 557 (1990). 

C. Itzykson and J.-B. Zuhei, Quantum Field Theory, pg 

276, McGraw-Hill, New York (1980). 

B. L. Altshuler, L. B. loffe, and A. J. Millis, Phys. Rev. B 
50, 14048 (1994); B. L. Altshuler, L. B. loffe, A. I. Larkin, 
and A. J. Millis, Phys. Rev. B 52, 4607 (1995). 

Y. B. Kim and A. J. Millis, Phys. Rev. B 67, 085102 
(2003). 



S. Sachdev and N. Read, Int. J. Mod. Phys. B 5, 219 
(1991). 

A. V. Chubukov, T. Senthil and S. Sachdev, Phys. Rev. 

Lett. 72, 2089 (1994); A. V. Chubukov, S. Sachdev, and 

T. Senthil, Nucl. Phys. B 426, 601 (1994). 

S. V. Isakov, T. Senthil, and Y. B. Kim, Phys. Rev. B 72, 

174417 (2005). 
^° S. Sachdev, Rev. Mod. Phys. 75, 913 (2003). 

P. Azaria, B. Delamotte, and T. Jolicoeur, Phys. Rev. Lett. 

64, 3175 (1990). 
^'^ D. P. DiVincenzo and E. J. Mele, Phys. Rev. B 29, 1685 

(1984). 



